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CN ■ Abstract 

We stochastically quantize the Born-Infeld field which can hardly be dealtwith by 
means of the standard canonical and/or path-integral quantization methods. 2 We 
set a hypothetical Langevin equation in order to quantize the Born-Infeld field, 
following the basic idea of stochastic quantization method. Numerically solving this 
nonlinear Langevin equation, we obtain a sort of "particle mass" associated with 
the gauge-invariant Born-Infeld field as a function of the so-called universal length. 
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O ■ 1 Introduction 



Many years ago Born and Infeld [1] presented a nonlinear electromagnetic 
field with a non-polynomial action including the so-called universal length. 
One of the most important characteristics of the Born-Infeld field is found in 
its static solution which has no divergence of static self-energy. Many physi- 
cists expected that this might be an example of divergence-free field theory. 
However, no one could succeed to quantize the field, by means of the standard 
canonical quantization method, because of the complicated nonlinearity. Even 
the path-integral quantization can hardly be applied to this field, because we 
cannot easily manipulate such a non-polynomial action. We have to invent a 
new quantization method if we want to quantize the Born-Infeld field. 

About ten years ago, Parisi and Wu [2] proposed a new quantization method, 
called stochastic quantization, by introducing a hypothetical stochastic pro- 
cess with respect to a new (fictitious) time, say t, other than ordinary time, 
say xq. The stochastic process is so designed as to yield quantum mechanics 
in the infinite t-limit (thermal equilibrium limit). This theory starts from a 
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hypothetical Langevin equation for the stochastic process, given by adding the 
fictitious-time derivative and the random source to the classical equation of 
motion. Remember that the stochastic quantization can be formulated only on 
the basis of classical field equation, without resorting to canonical formalism. 



2 Brief review of Born-Infeld field 



The ordinary electromagnetic field is described by the following Lagrangian 
density 

C = - 1 -F, V F^= 1 -(E"-U"), (1) 



where F^ v = d^A v — d u A^. The corresponding action is given by 

S = j d 4 xC . (2) 

Here we have followed the usual notation. Note that we keep the Minkowski 
metric. 

In the case of spherically symmetric static electric field, we can put 

E = -VA.(r), H = 0, (3) 

and obtain 

A = 6 - (4) 

for a point charge e. That this field becomes oo for r — > is a well-known fact. 
Let us introduce the action functional of the Born-Infeld field [1] : 



/ 



(fx 



(5) 



where 1/Vb is a sort of universal constant called universal length, and has the 
dimension of length in natural unit h — c — 1. We can easily see that Sb — * <S, 
that is, the Born-Infeld field will become the ordinary electromagnetic field as 
b tends to oo. 
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As is well known, all physical quantities can be written only in terms of the di- 
mension of length in natural unit h — c — 1. Many years ago Heisenberg antic- 
ipated that we could formulate a finite field theory, free from field-theoretical 
divergences, if we could bring a sort of universal length into physics in an 
appropriate way. Following his idea, Born and Infeld [1] proposed to use the 
above field given by action (5). Unfortunately, however, the Heisenberg's an- 
ticipation was not accomplished yet even now. 

For the spherically symmetric static field A (r) (3), the Born-Infeld action (5) 
yields the following equation 



d 
dr 



l-UiMr)) 2 , 



0, 



(6) 



whose solution is given by 



A = - f <% 2—^ 

r/ro 



(e/r) for r > r , 
1.8541 • (e/r ) for r -> , 



(7) 



where ro = y\e\/b. We surely realize that the Born-Infeld field has finite 
static self-energy. Of course, the self-energy goes back to the original infinity 
as b — > oo. 



We are not interested in static field but in wave field propagating to remote 
places. The Euler-Lagrange equation of Sb is generally written as 



8Sb 

5A v {x) 



■p\xv 



1 4- J-F 2 
1 + 262 r . 



X 



o , 



(8) 



or 



^F^ = (^ln X )i^, 



(9) 



where F 2 = F lul F fil ' and x 



21, 



The right-hand side of this equation is proportional to l/b 2 , and can be ex- 
panded in a power series of l/b 2 . Needless to say, its unperturbed one (for 
1/6 = 0) is nothing other than the free Maxwell equation 



d,F^ = , 



(10) 
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which allows us to use the wave gauge given by 
A = , V • A = . 



(11) 



Consequently, one may naively expect to have the perturbative theory based 
on (9) and (11). In this case, however, we can hardly develop this kind of the 
perturbative approach to the quantized Born-Infeld field, because the interac- 
tion part includes higher powers of derivative terms. This is the reason why we 
attempt to develop an unperturbative approach to the quantized Born-Infeld 
field by means of stochastic quantization [2] in the present paper. 



3 Stochastic quantization 

As is well-known, it is convenient to use the Euclidean action Se derived 
by the Wick rotation (xq — > — ixo) from the original Minkowski action, for 
the purpose of carrying out the Parisi-Wu stochastic quantization. In order 
to perform stochastic quantization of the Born-Infeld field A^, we have to 
introduce the additional dependence on fictitious-time t (other than ordinary- 
time xq) into the field quantities and then to set the basic Langevin equation 



for a hypothetical stochastic process of t) with respect to t [2]. Here 

rj^(x,t) is the Gaussian white-noise field subject to 



for statistical ensemble averages, where we have put h = 1 for simplicity. 

According to the prescription of stochastic quantization [2], we can derive the 
field-theoretical propagaters through the well-known formula 



where A^(x, t) as a function of is to be obtained by solving (12). 

Let us decompose A^ and rj^ into their longitudinal and transverse compo- 
nents, A L and A T , and rfc and rrf , given by 




(12) 



< ?7 M (:M) > = , 
< VtM (x, t)r) u (x', t') > = 25^5(x - x')5(t - t') 



(13) 
(14) 




(15) 
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Al = -d,d v A v (16) 

Al = (8^-^dpd v )A v - d,Al = 0, (17) 

and similar ones for 77's. Therefore, we can decompose the basic Langevin 
equation (12) as follows; 



d 



A^(x,t)=0 + ^(x,t) , 



dt 



SSe 

saj; 



\AT=AT(x,t)+Vfa,t) . 



(18) 
(19) 



The absence of drift force in (18) is an important reflection of the gauge 
invariance that Se does not depend on Af\ Consequently, the longitudinal 
component A^(x,t) makes a random walk around its initial value Aj^(x, 0) = 
^<9 M 0(x), 4>(x) being a scalar field. As was discussed in detail in the case of 
non-Abelian gauge field [3], we must introduce gauge parameter a by taking 
average of over random fluctuations around zero ((f>(x)(f>(x') = a5(x — x')). 
Thus we obtain 

D^(x,x') = ±dMa±- + 2t)5(x-x') . (20) 



Of course, we know that the longitudinal components never appear in gauge- 
invariant (physical) quantities. 

We should also notice that the transverse components and their propagaters 
are completely decoupled with the longitudinal ones in the present case. This is 
an important point quite different from the non-Abelian gauge field case. Thus 
we can safely discard A^, and use the wave gauge (11), even in the present 
case, for the purpose of deriving propagaters of the transverse components. 

Considering wave propagation along the z = rr 3 -axis, we put 



A = , A 3 = , 

Ax = Ai(x , x 3 ) , M = A 2 (x , x 3 ) , 



(21) 
(22) 



and 



7?o = 0, 773 = 0, (23) 
Vi = Vi(xo,x 3 ) , 7/2 = 7/2 (x ,x 3 ) . (24) 

Note that A^ and 7/ M have no longitudinal components. 
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In this case, the Euclid action becomes 
Se = J dx 



1 



(25) 



where F| = (<9 Ai) 2 + (d 3 A 2 ) 2 + (<9 A 2 ) 2 + (<9 3 Ai) 2 . The corresponding classical 
field equation is given by 



5Si 



SAi(x) 



= -d n 



d A t 



d 3 



d 3 Ai 



1,2). (26) 



Note that the dimensions of the field and b are different from the original ones: 
[Aj\ = [L°] and [b^ 1 ] = [L], [L] standing for the dimension of length. 

Based on the classical equation, we can set up the basic Langevin equation 
for stochastic quantization of the Born-Infeld field as follows; 



d_ 

dt 



Ai(x ,x 3 ,t) = d 



d A t 



+ d 3 



d 3 A t 



1 + ^1 



+ r)i(x ,x 3 ,t) , 



(i = l,2), (27) 
where the fluctuating source-field rji should have the statistical properties 

< rji(xo,x 3 ,t) > v =0 , (28) 

< 77 i (x ,x 3 ,i)77 J (xo,4,0 >r,= 26ij6(x - x' )S(x 3 - x' 3 )S(t - t') . (29) 

Consequently, we have to obtain Ai as a function (or functional) of rji, by 
solving (27), and then to calculate expectation values of physical quantities, 
by making use of (28) and (29). For example, the field-theoretical propagater 
of Ai is given by the formula 



AfAxo - x' , x 3 - x' 3 ) = lim [< Ai(x , x 3 , t)Aj(x' , x' 3 , t) > 



< Ai(x ,x 3 ,t) >< Aj(x' ,x' 3 ,t) >] . 



(30) 



For conventional fields, we can extract real information about the particle 
mass, M and M', associated with the field (or the first energy gap) from the 
asymptotic formulas 



A£(0,x 3 ) const. exp[-M|x 3 |] , 



(31) 
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Af t (x ,0) lxo -tT const. exp[-M'|x |] (32) 

(i: no summation), in which we can put Ai = Ai' for the Euclidean symmetry 
in space-time. Unfortunately in the Born-Infeld case, however, we have no 
reliable theory to justify the procedure (31) and/or (32) to give mass. Despite 
of this situation, we intend to follow the conventional approach to the "particle 
mass" associated with the (transverse) Born-Infeld field, based on (31) and/or 
(32). 



4 Numerical simulation and particle mass 



Needless to say, we know that it is very difficult to solve (27) analytically, so 
that we are inevitably enforced to deal with it by means of numerical simu- 
lation. For this purpose, we first discretize the Langevin equation (27) on an 
N x N lattice with spacings Axo and Ax 3 (along time and space directions, 
respectively). Denoting the Born-Infeld field on the (k,l)-th lattice point by 
Ai;k,i(f)i where i stands for the i-th component of the field and k, I for the or- 
dinary time and spatial position, then we write down the discretized Langevin 
equation as 

Aj. Al (t + At) -A rAl (t) 
At 



Gj-k+i,i(t) - Gj-k,i(t) Hj. A i + i(t) - Hj. A i(t) I 2 

Ax Ax 3 + V Ax Ax 3 At r ' kA } ' 

(33) 



where 



rr U \ ^(AiMt) ~ A i;k,l-l{t)) 

Hi;k,l{t) = 1 (35) 

with 



+ 



Al; k ,l{t) - A 1 . k _ l!l (t) 



AXn 



+ 



A 2 ;kM -^2;fc-l,l(*)' 



Ax 



~A 1;k ,l(t) - Al; fc|i -l(t) 

~~ Ax, 



1 2 



+ 



A 2 ;k,l( t ) - A 2;k,l-l(t) 

Ar A 



1 2" 



(36) 
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for drift terms, and 



< N i;kj i(t) > N = 0, < Ni. !k> i(t)N i . k > jl >(t') > N = Sij5 kk >5u'5 tt ' , (37) 
= 1,2), (l,k = l,...,N) 

for noise terms. 

Here let us introduce a scale unit a, which has dimension of length, and put 
relevant quantities in the following way: 

Ax = Ax a, Ax 3 = Ax 3 a, At = Ata 2 , b = ba -1 , 
Ai;k,i(t) = Ai; k ,i(t), Gi; k) i(t) = G i;k; i(t)a 1 , 
H Mk ,i(t) = H rAl (t)a-\ N Mk:l (t) = Ni- k ,i(t), 

M=Ma~\ (38) 

Note that all quantities with tilde are dimensionless. Thus, the equations (33) 
and (37) are rewritten as 

Aj. kt i{i+ At) -Ai^ijt) 
At 



_ Gj. k+1: i(t) - Gj. A i{t) Hj. Kl+l {t) - Hj. Kl {t) / 2 ~ ~ 
Ax + Ax 3 + \j Ax Ax 3 At r ' kA ) ' 

(39) 

< Ni-k,i(t) >ff= 0, < Ni- k ,i(t)Ni. k/ , v (t') >x= SijSwSwStt, . (40) 

Gi ;kt i and Hi ]kt i include the dimensionless _F| in the same way as Gi- k j and 
H i]kt i depend on F|. Note that only F| contains the Born-Infeld parameter b. 

In the conventional field theory, this scale unit a can be determined by making 
use of the renormalization group theory. In the present case, however, we have 
no reliable theory to determine the scale unit, and then we shall inevitably 
calculate all quantities (in particular, the "particle mass") on an arbitrary 
scale. Only for the sake of simplicity, let us put a = 1, in order to go on 
our procedure. For a while from now on, we suppress those tilders which are 
put on the quantities. Therefore, the "particle mass", M., will be given as a 
dimensionless quantity in this scheme. 

In order to solve numerically the above equation and obtain the field-theoretical 
propagater Ajj(x , x 3 ) on the above lattice, we should introduce the periodic 
boundary condition given by 

Ai ]k ,i+2i c = A i;kt i , Ai. k+2 k c ,i — A i;k; i , (41) 



8 



Fig. 1. Field-theoretical propagaters for several values of b . 

where l c , k c and 2Z C , 2k c stand for the lattice center and the period, respectively. 
Practically, we have used the Langevin-source method (for example, see [2]), in 
which we have performed 5.4 x 10 6 iterations for a lattice of 20 x 20 sites with 
At = 0.01 and l c = k c = 10 (to realize thermal equilibrium), and then use the 
subsequent 2.0 x 10 5 iterations to calculate the field-theoretical propagater. 

Figure 1 shows our numerical results of the field-theoretical propagaters An(0, X3) 
for b^ 1 = 10,20,30,40 and 50, together with curves given by the asymptotic 
formula 

AJX, 0) = AJ0, X) = c coshM \X -*c\ ( , . no summa ti on ) ? (42) 

coshAl|x c | 



where x c stands for the center of lattice. Also we numerically have got similar 
field-theoretical propagaters for other directions and/or components. Equa- 
tion (42) is the substitute of (31) and/or (32) under the boundary condition 
(41). We have estimated the "particle mass", Ai, associated with the Born- 
Infeld field, by making use of x 2 -fitting based on (42). As shown in Table 
1, our results are the following: M = 0.0415,0.0835,0.1260,0.1693,0.2162, 
correspondingly to b' 1 = 10,20,30,40,50, where x 2 is 4.69 x 10~ 4 . We have 
estimated the statistical fluctuations for the fictitious time as accuracies in 
Table 1. Here we have put C = Ajj(0,0) =< Af > (i: no summation). Note 
that C is independent of i due to the space-time uniformity, and that C is 
gauge-invariant . 

Rigorously speaking from the point of view of (38), we can only assert that the 
above M. is proportional to the "particle mass" associated with the (trans- 
verse) Born-Infeld field. We should repeat that we have no renormalization 
group theory to give the scaling formula in the case of Born-Infeld field. Re- 
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Table 1 



The "particle mass", A4, associated with the Born-Infeld field 



b- 1 


10 


20 


30 


40 


50 


M 

Accuracy 


0.0415 
±0.0062 


0.0835 
±0.0021 


0.1260 
±0.0017 


0.1693 
±0.0017 


0.2162 
±0.0015 



Fig. 2. "Particle mass", M., as a function of b 



Fig. 3. "Particle mass", M.', as a function of b 1 . 

member that the problem is still open to questions. In this paper, however, 
we are talking about the "particle mass" by M. which is given by (42). 

Figures 2 and 3 plot the "particle mass", given by M. and M, as a function of 
b^ 1 . Observe that M = M'. It seems that the "particle mass" is proportional 
to but unfortunately, we do not know what kind of physical implications 
this fact suggests. Another important point should be that the "particle mass" 
seems vanishing, in the case of b^ 1 = 0, as expected from the fact that the 
Born-Infeld field must go back to the free Maxwell field in this limit. All results 
are presented in Table I. 
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Table 1 also tells us that the (dimensionless) "particle mass" on this scale 
(with a = 1) distributes over a very small region. 

Table 1 or Figures 2 and 3 can be fitted well by a single formula given by 

M = 77 , 7 = 0.00426 . (43) 
b 



This equation is rewritten in terms of M. and b 1 having dimension as 

M = (^)j, 7 = 0.00426. (44) 
cr o 



Here let us try to choose a = b 1 as the scale unit, then we obtain 

M — jb, 7 = 0.00426 , (45) 

or 

M = 7, 7 = 0.00426 , (46) 

in the dimensionless expression, because b = 1 in this case. This implies that 
the constant 7 is nothing other than the "particle mass" , being independent 
of the universal length, on the scale adjusted by a = b~ x . 

Finally, let us examine whether our "particle mass" M. can be regarded as a 
sort of particle mass in the sense of conventional field theory. For this purpose, 
we should numerically compute Fourier transform h.a{k 2 ) given by 

N r dr N r dr 

Adk 2 ) =Jy=J ^e-^°-^3 A ^ 3?Xo) > (47) 
v v 



iV standing for lattice size. If M. meant a sort of particle mass in this sense, 
we could hardly observe so sharp .M-dependence of hu(k 2 ) as a function of 
k 2 , for v / ^> M = 0.0415,0.0835,0.1260,0.1693,0.2162. Figure 4 shows our 
numerical results (see [4] for technical details), in which all curves for various 
fe's are normalized to Ajj(0) = 1. We can clearly observe in Fig. 4 that Aa(k 2 ) is 
almost independent of k 2 , except in its height. That is to say, our anticipation 
seems justified. For comparison, we put the dashed curve representing a free- 
propagater (being a substitute of massless free Maxwell field), e 2 /(k 2 + e 2 ), 
with a very small mass e = 0.001 M. 

In order to reconfirm this fact, we present Fig. 5, (for ^u(k 2 ) versus k 2 ) 
stressing that all curves overlap each others for larger k 2 . Note that long 
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Fig. 4. A n (fc 2 ) corresponding to masses obtained in Table 1 



Fig. 5. -^pAii(k 2 ) versus k 2 

distance behavior of the propagater in configuration space gives A4, while 
short distance behavior determines Fourier transforms for larger k 2 . 

Moreover, we compare a special Aa(k 2 ) with the corresponding Feynman prop- 
agater M 2 j{k 2 + M 2 ) for M. = 0.126 in Fig. 6, as an example. In this figure 
one could observe that the difference between them would represent a possible 
(damping) effect due to the non-linearity of the Born-Infeld field. 
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Fig. 6. Comparison of An(/c 2 ) with M 2 /(k 2 + M 2 ) for M = 0.126. 
5 Conclusion 

Summarizing, we have stochastically quantized the Born-Infeld field, charac- 
terized by the so-called universal length, which cannot be dealt with by means 
of the conventional quantization methods. Even though we can hardly justify 
the whole procedure theoretically, we have derived the "particle mass" asso- 
ciated with the (transverse) Born-Infeld field, as a function of the universal 
length, through the conventional formulas to give them. 

It would be interesting to observe that we have derived the "particle mass" 
from a perfectly gauge-invariant field theory. Of course, we can guess that the 
"particle mass" is produced by introducing the universal length b^ 1 having 
the dimension of length. 

The authors are indebted to Drs. I. Ohba, S. Tanaka, Y. Yamanaka, K. Okano 
and B. Zheng for many discussions and suggestions. 
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